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Abstract. We prove that the unique least-perimeter way of partitioning the unit 2-dimen- 
sional disk into three regions of prescribed areas is by means of the standard graph described 
in Figure 



Introduction 

Partitioning problems in the Calculus of Variations have multiple applications in physical 
sciences. They can model multitude of natural phenomena such as the shape of a cellular 
tissue, the interface of separation between fluids, and many others, as described in the treatise 
by D'Arcy Thompson |Thj . 

In this work we consider the isoperimetric problem of partitioning a planar disk into three 
regions of given areas with the least possible perimeter, and we prove that the standard 
configuration in Figure Q consisting of three circular arcs or segments meeting orthogonally 
the boundary of the disk, and meeting in threes at 120 degrees in an interior vertex, is the 
only solution to this problem. 




Figure 1. The least-perimeter partition of the disk into three given areas 

In addition to the above conditions, the solution must satisfy a certain balancing condition 
on the geodesic curvatures of the circles. This condition will be stated precisely in the next 
section. 

Existence and regularity of solutions for this problem are guaranteed by the results of 
F. Morgan M2], who showed that the minimizer, in the interior of the disk, is composed 
of smooth curves of constant geodesic curvature meeting in threes at 120 degree angles. 
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Boundary regularity also follows from M2 although it is not explicitly stated in his work. 
Existence and regularity in higher dimension were studied by F. Almgren |Almj . 

The least-perimeter way of partitioning a disk D into two regions of given areas is by 
means of an arc of circle or segment that meets orthogonally dD. From the existence and 
regularity results in next section it follows that there is a solution, which is a smooth, possibly 
nonconnected, embedded curve with constant geodesic curvature that meets dD orthogonally. 
Such a curve must be connected, since otherwise we could rotate one component with respect 
to the center of the disk until it touches a second one, thus producing a non-allowed singularity. 
On the other hand, as the curve has constant geodesic curvature, it must be part of a circle 
or of a line. 




Figure 2. The least-perimeter partition of the disk into two given areas 

The isoperimetric problem consisting of enclosing n given areas in the disk or in the plane 
with the least possible perimeter has a strong complexity which is derived not from the 
geometry of the individual components of the solution (they can be described in terms of 
circles or lines) but from their large number. 

The planar double bubble conjecture was proved by J. Foisy et al. FABH , who showed in 
1993 that the standard planar double bubble uniquely minimizes perimeter in M 2 . Assuming 
that the studied regions are connected, C. Cox et al. |CHH| proved in 1994 that the standard 
planar triple bubble uniquely minimizes perimeter in the plane for any three given areas. 
R. P. Devereaux [O] studied in 1998 the planar triple bubble conjecture under the hypothesis 
that all the regions have the same pressure. W. Wichiramala finally proved the planar triple 
bubble conjecture in 2002 in his Ph. D. Thesis |Wj . J. Masters |M| proved in 1996 the double 
bubble conjecture in S 2 . Interesting preliminary work was carried out by Bleicher |Blj . |B2j . 
|B3j . Concerning boundary problems, G. Hruska et al. |HLPSj have obtained some results 
for planar bubbles in corners. Also results on tori and cones have been obtained in C HLLj . 

m 

In higher dimensions, J. Hass and R. Schlafly |HSj proved the double bubble conjecture in 
R 3 for equal volumes. The general conjecture was settled by M. Hutchings et al. HMRR . 
For higher dimensional Euclidean spaces Reichardt et al. RHLS have obtained a proof of the 
double bubble conjecture in R 4 and partial results in higher dimensional Euclidean spaces. In 
the three-dimensional torus, M. Carrion et al. |CCWB] have provided numerical evidence for 
a double bubble conjecture with ten types of solutions. In the three-dimensional sphere and 
the three-dimensional hyperbolic space, A. Cotton and D. Freeman |CF| have also obtained 
partial results on the conjecture that the standard double bubble in these spaces uniquely 
minimizes perimeter. 
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Planar bubbles are also of great interest to physicists. Interesting articles focusing on 
physical aspects of the problem are |G.I.IF| and CGFM j. 

The most interesting mathematical open question for these problems is to show that the 
minimizing configurations must have connected regions, either in M? or in the disk. In ad- 
dition, in the planar problem, one should also be able to prove that the exterior region is 
connected, i.e., that there are no empty chambers. 

We have organized this paper in several sections. In Section ^ we give precise definitions, 
compute the first and second variations of length for graphs, recall existence and regularity 
results for the problem of minimizing perimeter while partitioning the disk into given areas, 
and state some properties which minimizing graphs must satisfy. In Section [2] we obtain a 
bound on the number of components of the largest pressure region determined by a graph 
which minimizes perimeter up to second order. We conclude that a minimizing configuration 
must have one often possible types, described in FigureHJ In Section|2]we prove the necessary 
results to discard the possibilities obtained in Section which allow us to prove our Main 
Theorem in Section^ In a final section, we indicate further lines of research and give several 
conjectures. 

All the pictures in this paper have been made by using Surface Evolver, a software developed 
by Ken Brakke (http://www.susqu.edu/facstaff/b/brakke/). 



1. Preliminaries 

1.1. Notation. Let D C M? be the closed unit disk in centered at the origin. An admissible 
graph C C D consists of vertices and curves so that at every interior vertex (that is, a vertex 
in the interior of D) three curves of C meet and at every boundary vertex (a vertex in dD) 
just one curve of C meets dD. We shall also assume that C induces a decomposition of the 
open unit disk into n regions Ri, 1 ^ i ^ n, possibly nonconnected. An m-component is a 
connected component of a region with m edges. 

If R. L and Rj are adjacent regions, we will denote by Cij C C the (not necessarily connected) 
curve separating them. Let = {j ^ i;Rj touches Ri}. With this notation 

dRi n int(iJ) = |J dj. 

We shall denote by Nij the normal vector to the curve Cij pointing into the region Ri, and 
by hij the geodesic curvature of the curve Cy with respect to the normal . 

A standard graph consists in three circular arcs or lines segments meeting at an interior 
vertex at 120 degree angles, reaching orthogonally dD, and so that the sum of the geodesic 
curvatures is zero. 

Given n positive numbers a±, . . . ,a n such that Ya=1 a i = 7r ' the isoperimetric profile is the 
function I{a\, . . . , a n ) defined as the infimum of the lengths of all admissible graphs separating 
regions in the disk of areas a\, . . . , a n . 

We will say that an admissible graph C is minimizing for prescribed areas a\,...,a n if 
I(a\, . . . , a n ) is attained by C. 
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1.2. Variational formulae. Given an admissible graph C C D, we will consider smooth 
one-parameter variations (ft : C — > D for £ small, which satisfy <pt(dD) C 3D. We will denote 
by X = dip t I dt\ t= Q the associated infinitesimal vector field, which is smooth on every curve 
Cij. Note that X{p) is tangent to dD for each p in dD. Let = X ■ Nij be the normal 
component of X on Cy. 

Given such a variation, it is easy to check that the derivative of the area Ai of R4 at t = 
is given by 



*=« jG/(i) C ' 



For the derivative of length for such a variation we have 

Proposition 1.1 (First variation of length HMRR, Lemma 3.1]). Consider an admissible 
graph C C D, and a smooth variation (ft : C — > D with associated vector field X . Then the 
first derivative of the length of <ft(C) at t = is given by 



< L2 » 



i=0 ie{i,...,n} pedd. 



where Vij{p) is the inner conormal to Cij in p . 

We will say that an admissible graph is stationary if (|1.2|) vanishes for any area-preserving 
variation. From Proposition II. II it is easy to prove the following 

Proposition 1.2. Let C <Z D be a stationary graph. Then the following conditions are 
satisfied 

(i) The geodesic curvature hij is constant on C^. 

(ii) The edges of C meet in threes at 120-degree angles in interior vertices. 

(hi) The balancing condition: three edges C^, Cjk, C^i meeting in an interior vertex satisfy 

(1.3) h^ + hj k + hki = 0. 

(iv) The edges of C meet dD orthogonally at boundary vertices. 

Condition (ii) implies that, in some interior vertex where the three curves Cij, Cjk, Cki 
meet, the normals add up to zero, i. e., Nij + iV^ + N^i = 0. This implies that the normal 
components of the vector field X must satisfy 

(1.4) Uij + Uj k + u ki = 0. 

Given a stationary graph C, and a function u : \J i ■ Cij — > R, with mj = w|c y , satisfying 
condition (|1.4|) on every interior vertex, it is always possible to find a vector field X on 
C, so that Uij = X ■ Nj and X is tangent to dD in each boundary vertex. Associated 
to X one can also find a one-parameter variation tpt : C — > I?, for £ small enough, so that 
<Pt(p) = ex P p (^(p)) for any p out of an arbitrarily small neighbourhood of dD. The argument 
is as follows: fix an arbitrary neighbourhood U of dD that does not contain interior vertices 
of C. Modify X so that it is normal to C in U. Let v be the inner normal to dD. Extend 
it to U so that it is tangent to the edges of C. Also extend X to a vector field on U D D 
by means of the exponential mapping. Let A be a smooth function equal to 1 near dD with 
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support in U. Consider the vector field Y = X — (X ■ (Xu)) \v and the local one-parameter 
group -0* generated by Y. Since v is tangent to C and X is normal to C in U, we have that 
Y = X on C. Moreover, for p e dD, the vector Y{p) is tangent to dD. Hence the deformation 
ipt(C H U) has initial velocity vector field X and keeps C inside the disk. The variation ipt(C) 
has the further property that coincides with exp p (tX (p)) in U out of the support of A. Now 
we simply define (ft(p) equal to exp p (tX(p)) out of the support of A, and equal to 4>t(p) in U. 

The balancing condition (|1.3|) allows us to define a pressure pi on every region Ri, starting 
from a given region, so that 

(1.5) h ij =p i -p j . 

These pressures are determined up to an additive constant. The first variation formula of 
length can be rewritten in terms of pressures in the following way: if C is a stationary graph, 
then the first variation of length for an arbitrary variation is given by: 



(1.6) 



dL 
HE 



EdA.i 

i=i 



Observe that the indetermination of the pressures up to some additive constant does not affect 
the above formula since Ya=i dAi/dt = for any variation of the regions Ri, as J27=i Ai(t) = 
area(-D) along the variation. 

Let us prove now the second variation formula of length 

Proposition 1.3 (Second variation of length). Let C be a stationary graph and let {<pt} be a 
variation with associated vector field X preserving areas up to second order. Then the second 
derivative of length at t = is given by 



(1.7) 



E 

i=l,...,n 



hf j u ij )u ij + 

pGint(D) 



dUi 



(p) 



+ 



E 

pedD 



u 2 , 



where qij(p) = (hk% + hkj)(p)/V3, and R^ is the third region touching the vertex p. 



(P) 



Proof. Differentiating the integral terms in equation (|1.2[) . we get 

( f hijUij] = j (u"j + hf jUij ) Uij + hij ^ 
dt t=0 \J Cij J Jc l3 dt t=0 \JC i3 

but since the variation preserves areas up to second order, it follows that 



a 



d ( f \ n 
E h vT f ( L Ui i) = 2 E^ 

l,...,nl t=0 V Jc ^ ' j=i 



ie{l,...,n} 

ie/(i) 



d 2 A { 



dt 2 



0. 



t=o 



Differentiating now the second term in equation (|1.2|) . we get 
d 



dt 



(X ■ Uij) = (D X X ■ Uij) + mj h^ (X ■ u^) + Ui 



duj 
' dui 
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For p £ int(-D), the first term vanishes since v, L j + Vjf. + = 0, and after some calculations as 
in HMRR , the second one can be seen as —qijufj, where = (hki + hkj)/V3- For p G dD, 
since the configuration is stationary, the edges meet dD orthogonally, so that DxX(p) ■ Vij{p) 
equals ufj times the geodesic curvature of dD, and (X ■ Vij){p) = 0. □ 

The condition that the variation must preserve area up to second order is not really needed 
as we can show in the next Lemma 

Lemma 1.4. Let C C D be a stationary graph. Given smooth functions Uij : Cij —* IR such 
that (jl.ljl and (|1.4j) are satisfied (a variation that preserves area up to first order is given), 
there is a variation {(ft} of C which leaves constant the area of the regions enclosed by (ft(C) 
and such that the normal components of the initial velocity vector field X are the functions Uij . 

Proof. Let X be a vector field on C, smooth over each curve Cij, such that X ■ Ny = Uij. Let 
ip t : C — > D be a one-parameter variation of C associated to X such that iptip) = exp (tX(p)) 
out of a small neighbourhood U of dD which does not contain interior vertices of C. 

We label the regions Ri so that Ri touches Ri+i for i = 1, . . ., n — 1. Choose positive 
functions Vi with support in the interior of C^j+i) and out of U. The variation induced by the 
vector field ViN^ i+1 ^ decreases the area of Ri, increases the area of Ri+i and leaves constant 
the area of the remaining regions. 

Consider the variation equal to 

(t,3i,.. . ,s„_i) i — ► exp p (tX(p) + ^2 s i v i N i(i+i)(p)) , in Cn (£> - U), 

i=l 

and equal to ipt(p) for p £ C HU. Consider the function (A±, . . . , ^4 n _i) of (t, s%, . . . , s n _i), 
given by the areas of the deformation of the regions R\,. . ., R n -i- The Jacobian 

d(A 1 ,...,A n -i) 
d(si,... ,s„_i) 

is lower triangular, with non-vanishing entries in the principal diagonal, so that the matrix is 
regular. The Implicit Function Theorem allows us to find smooth functions s\(t),. . ., s n _i(t) 
such that Ai(t, si(t), . . . , s n _i(i)) is constant for all i. 

The initial velocity vector field of such a variation is equal to X on C fl U, and to X + 
J27=i s i(0) ViNtfi+i) on C fl (D — U). As s-(0) = since ipt preserves areas up to first order, 
we conclude that X is the initial velocity vector field. □ 

Remark 1.5. A variation of a stationary graph C by stationary graphs preserves the angles 
between edges at interior vertices and the orthogonality condition at boundary vertices. Given 
a variation preserving the area of all the regions up to first order, we can modify it by 
Lemma 1 1.41 so that the areas enclosed are constant along the deformation. From the second 
variation formula we get that the second derivative of length is given by 

d 2 L ^ dPa dA a 

where a labels the components of the stationary graph (regions can be disconnected), and 
dp a /dt is the derivative of the pressure of the component a with respect to the considered 
variation. Take into account that the quantity u'(, + hf-Uij, the derivative of the geodesic 
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curvature h™, only depends on u^j, the normal component of the variational vector field X, 
and that the modification needed in Lemma fl.4l to preserve areas only modifies the acceleration 
of the variation. The angle-preserving condition depends only on the initial velocity vector 
field. 



In general, if the areas are not preserved up to second order, the second derivative of length, 
for a deformation of a stationary graph by stationary graphs, is given by 

d 2 L ^ dp a dA a d 2 A a 

n? = ^~dr^jr +Pa ~dW : 

a 

which can also be obtained by differentiating equation (|1.6jl . 

Remark 1.6. For a variation such that the angles between the edges are preserved, we have 
Dxi^ij + Vjk + Vki) = (since i/y + Ujf. + = for all t), so the boundary term in the second 
variation formula vanishes. 



1.3. Admissible functions and the index form. Let C be a stationary graph. We say 
a function u : UijCij — > IR. is admissible if the restrictions Uij = u\dj he in the Sobolev 
space W 1,2 (Cij), and verify that at any interior vertex p, Uij(p) + %fe(p) + Uki(p) = 0. These 
functions correspond to variations of C which have as normal components of the associated 
vector field X the functions Uj,-. These variations will preserve areas if, for each i, 



E L ^ = - 

An admissible function u is said to be a Jacobi function if the associated variation preserves 
the geodesic curvatures of each edge Cjj and the angles in each vertex. The fact that the 
geodesic curvatures are preserved means that the restrictions Ujj to Cjj verify 

u ij + tfjUij = 0. 

It is clear that the normal component of the Killing vector field generated by the rotations 
about the origin gives a Jacobi function. 

From equation (jl.7|) . we define the index form, that is, the induced bilinear form defined 
on the space of admissible functions, by 

(1.8) Q(u,v) = --l J {u"j + h 2 jUij ) 

i=l,...,n *i 

ie/(i) 

peint(D) pedD 

where qij are the functions defined in Proposition 11.31 

We will say a stationary graph C is stable if Q(u,u) ^ for any admissible function u 
whose associated variation preserves areas, and unstable if it is not stable. It is clear that a 
minimizing configuration must be stable. 
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1.4. Existence and Regularity. From the results of F. Morgan |M2j . we obtain the fo- 
llowing 

Theorem 1.7 (Existence and Regularity |M21 Th. 2.3]). Let D Cl 2 be a closed disk, and let 
a\, . . . , a n be n given areas such that Y^i=i a i = area(D). Then there exists a graph separating 
D into n regions of areas a±, . . . ,a n . Moreover such a graph consists of constant geodesic 
curvature curves meeting in threes in the interior of D at 120 degree angles, satisfying the 
balancing condition (jl.3j) for the geodesic curvatures, and meeting dD, one at a time, in an 
orthogonal way. 

Proof. From the results in |M2j one gets the existence of a solution and the regularity in the 
interior of the disk with just triple points as possible singularities. One also gets that there 
is a finite number of components (and hence of curves) in the minimizing configuration. For 
the boundary regularity, we only need to prove that at every point of dD, at most one curve 
of the minimizing configuration arrives, at 90 degrees. 

If one or several curves meet dD at p and at least one of them is not orthogonal to dD, then 
the first variation formula implies that the graph is not stationary. Suppose now that several 
curves meet orthogonally dD at p. We order them counter-clockwise and we consider the first 
one, C, which is the common boundary of components Qi and fL, with Qi a boundary one. 
Make a small deformation in the interior of C which implies a loss of area 5 for fij. In order 
to preserve the areas, it is possible to choose a point q near p in C, join q to the second curve 
C , which is in the boundary of Qj, and eliminate the part of C between p and q. It can be 
checked that this new configuration, for 5 small enough, reduces perimeter. Then we get the 
desired regularity in the boundary of D. □ 

1.5. Some properties of minimizing graphs. We now give and recall some results on 
minimizing graphs that will be used to prove our main theorem 

Lemma 1.8. Given n positive numbers a\, . . . , a n such that J27=i a i = 7r > we have 



Moreover, equality is never achieved for n ^ 4. If equality holds in the case n = 3 then the 
standard graph consisting of three line segments dividing the disk into three regions of equal 
areas is minimizing. 

Proof. We can divide the disk into regions of given areas a\, ... ,a n by using appropriate n 
radii. This gives Q1.9JI . For n ^ 4, this configuration has a prohibited singularity at the origin, 
so that it cannot be minimizing. If equality holds in Q1.9JI for n = 3, the configuration must be 
stationary, so that the three radii meet in 120 degrees, and the configuration is the standard 
one for equal areas. □ 

Lemma 1.9. A minimizing graph must be connected. 

Proof. On a nonconnected graph, we can rotate one of the components until it touches another 
one creating an irregular meeting, so the graph cannot be minimizing. □ 

Remark 1.10. Let C C D be a minimizing graph, and f2 a connected component of D — C. 
Lemma ll . 91 implies that d£l D dD has to be connected. 

Lemma 1.11 ( FABH, Lemma 2.4]). On a minimizing graph, there are no 2-components. 




I(a 1 ,.. . ,a n ) ^ n. 
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2. A BOUND ON THE NUMBER OF COMPONENTS OF THE LARGEST PRESSURE REGION 



Lemma 2.1. Let C be a stable graph separating D into n regions. Then the region of largest 
pressure has at most n — 1 nonhexagonal components. 

Proof. Assume R\ is the region of largest pressure and suppose it has at least n nonhexagonal 
components, Q\, . . . ,Q n . For each i, consider the variation given by Ui = 1 on 9f2,, extended 
by zero to the whole graph. If 0^ is a boundary component then 

Q( Ui , Ui ) = - Y, { I h lj + E -mj(p)+ E 1 }<°> 

peint(D) p€dD 

since, for p in Cy n C\k fl dfli, we have 

t \ , / \ ^fei + hkj + hji + hjk f v /j-fci + /iji / s / n 
3ii(p) + 9ifc(p) = —f^ 1 — (?) = ^ ^ °' 

as i?i has the largest pressure. 

If Vti is an interior component then Q(ui,Ui) can be computed as above except that the 
last summand does not appear. So we get 

Q(ui,Ui) < 0, 

and equality holds if and only if fij is bounded by segments. It is easy to obtain, from Gauss- 
Bonnet Theorem, that Qj has to be an hexagon. In the case of three regions, this only occurs 
if the three pressures are equal. 

Hence, in our case we can find some nontrivial linear combination u of Uj, such that the 
induced variation preserves areas up to first order and Q(u,u) < 0. □ 

Lemma 2.2. Let C C D be a minimizing graph separating D into three regions. Then C is 
one of the graphs in Figure ^| 

Proof. Suppose first that all the pressures are equal. If all the components touch the boundary 
of D then C is standard. If there is an interior component, then it is hexagonal. It is easy 
to see that the edges leaving the vertices of the hexagon meet dD (otherwise we could find 
two different parallel rays meeting orthogonally dD). This implies that the graph is like in 
Figure This graph has two regions with three nonhexagonal convex components, and so it 
is unstable by Lemma 12. II 




Figure 3. A graph with an hexagonal interior component 




Figure 4. The ten possible configurations for minimizing graphs 
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Assume now that p\ ^ P2 ^ P3, with p\ > p%. Then R± cannot have hexagonal components 
and so has at most two components, by Lemma 12.11 An interior component of R\ has an 
even number of edges in its boundary and cannot be a 2-component by Lemma 11.111 So 
any interior component of R\ is a quadrilateral. A boundary component of R\ will have by 
Gauss-Bonnet Theorem three or four edges in its boundary. 

Suppose R\ is connected. If it touches dD and has three edges, we have the standard 
configuration [(10)] and if it has four edges, we have configuration (1) If R\ is interior , it 
will have only four edges, corresponding to configuration |(3)| 



Suppose now R\ has two connected components, A and B. They can be interior or boun- 
dary components. We study each case. 

If both are boundary components, as before, they can have three or four edges. If A and 



B have three edges, by the connectedness of C, the only possibility is configuration (1) If A 



has four edges and B has three edges, the only possibility is configuration [(7)1 and if A and 



B have four edges, we will have configurations (2) and|(3) 



If A is an exterior component and B is an interior component, as above, A will have three 
or four edges, and B will have four. In the first case, we will get configuration |(4)[ and in 
the second one, configurations [(5)| and [(8)| 

If both components are interior, they will have four edges. A component of the region 
with the smallest pressure cannot be interior by Gauss-Bonnet Theorem. Then the only 
possibilities are configurations [(6)| and [(9)| □ 



3. Unstable and non-minimizing configurations 

Lemma 3.1. Let C12 C D be a circle or segment meeting dD orthogonally between two 
regions R[, R' 2 with associated pressures p\ and P2 {the geodesic curvature of C\2 w. r. t. the 
normal pointing into R[ equals p\ — P2Y 

Then, given v £ C12, there exist unique curves C23, C31 with constant geodesic curvature 
yielding a standard graph. Moreover, the pressure p% of the third region is a monotone function 
of v taking on all values between —00 and +00, and the regions R% are nested. 

Proof. Fix v G C\2- Let q be a point in C12 fl dD. Consider the Mobius transformation 

i(z + q) 



m 

q — z, 

which takes the disk D to the upper half-plane and sends q to infinity. Then /(C12) is a 
straight line L. Assume that there are C23, C31 curves with constant geodesic curvature 
meeting C\2 at v so that C\2 U C23 U C31 is standard. Since / is conformal, the sum of 
the geodesic curvatures of f(Cij) is zero. Moreover /(C23), /(C31) intersect the real axis 
orthogonally. It is not difficult to see that /(C23), /(C31) are circles centered at the real axis, 
with the same radius by the balancing condition (|1 .3|) . Thus /(C23), /(C31) are unique and 
so are C23 and C31. It is clear that the regions determined by /(C23) and /(C31) are nested, 
which implies that i?| are nested. 
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Let d be the distance from f(v) to the real axis. By applying the inverse of /, it is possible 
to compute the geodesic curvatures /131 and I132 in terms of d, obtaining 



Hence p% = p\ + h^x, which decreases from +00 to —00 when d moves from to +00. □ 

Proposition 3.2. Given three pressures, there is a standard graph separating the disk into 
three regions with the given pressures. Such a graph is unique up to a rigid motion of the disk. 

Proof. Given two pressures p\, P2, there is a circle or segment Cyi separating D into two 
regions with constant geodesic curvature p\ — P2 meeting dD orthogonally. C12 is unique up 
to a rigid motion of the disk. Using Lemma 13.11 we can find a vertex v € C12 and unique 
curves C23, C31 so that C12 U C23 U C31 is a standard graph separating the disk into three 
regions with the given pressures. Uniqueness follows from the construction. □ 

Proposition 3.3. Let C\, C2, C3 be circles or lines meeting at 120 degree angles at some inte- 
rior point of D, satisfying the balancing condition (jl.3|) . If C\ and C2 meet dD orthogonally, 
then so it does C3. 

Proof. Let be the region enclosed by C\, C2 and dD. Apply the disk onto the upper half- 
plane by means of a Mobius map / sending p E dD — dQ to infinity. It is enough to show that 
/(C3) meets orthogonally the real axis. As in the proof of Lemma l3~Tl f(C±) and /(C2) meet 
themselves at 120 degrees and the real axis at 90 degrees. As / is a Mobius transformation, 
the balancing condition (|1.3f) is preserved. In case /(C3) is a line, it meets the real axis 
orthogonally. If /(C3) is a circle, then f{C\) U /(C2) U /(C3) is a standard planar double 
bubble, for which the centers of the circles are known to be aligned. The centers of f{C\) 
and /(C2) lie in the real axis, and hence also the center of /(C3). So we conclude that /(C3) 
meets the real axis orthogonally. □ 

The proof of Lemma 13.11 establishes the existence of a deformation of one of the regions 
along one of the edges. More precisely we have 

Proposition 3.4. Given a stationary graph C with a boundary 3-component £1, there exists 
a variation of C that 

(i) strictly increases the pressure of £1 while keeping the other pressures unchanged, and 

(ii) strictly decreases the area ofVL, and 

(iii) leaves invariant the edges of C not lying in dtt . 

Proposition 3.5. Let C be a stationary graph in which a region has two boundary 3-compo- 
nents. Then C is unstable. 

Proof. Let Ox, ^2 be boundary 3-components of the same region. On each Qj, consider the 
variation given by Proposition 13.41 The normal components U{ of the associated variational 
vector fields have disjoint supports and satisfy Q(ui,Ui) < 0. By Prop ositions 13 . 21 and 13 . 31 Q\ 



(3.1) 




(3.2) 
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and ^2 are congruent so that u = U\ — U2 satisfies the mean value conditions Hence 
the graph is unstable. □ 

Theorem 3.6. Given three areas a\, 02, 03 such that 01 + 02 + 03 = area(D) ; there is a unique 
standard graph, up to rigid motions of the disk, separating D into three regions of areas a{. 

Proof. Consider two standard graphs: C = C12 U C23 U C31, and C = C 12 U C 23 U C 31 so that 
Oj = area(i?j) = area(i?-), i = 1, 2, 3. 

In case h\ 2 = h' 12 we can a PPby a rigid motion of the disk to C until C\i and C[ 2 coincide 
near dD. As the areas of the enclosed regions are equal, Lemma 13 . 1 1 implies that C = C . 

Assume that h\2 > h' l2 - By Lemma 13. II we can continuously decrease the pressure p\ (while 
keeping constant P2 and p$) until we get another standard graph C" with h'{ 2 = hyi- For 
this new graph C" we get area(i?'/) > a\, area(i?-') < en, for i = 2, 3. Now we can apply 
an isometry of the disk to C" so that C'{ 2 and C[ 2 coincide near dD. As area(i?'/) > ax, 
Lemma 13 . 1 1 implies that area(i? 2 ') > a 2, which gives us a contradiction. □ 

Given two edges in a graph C, we will say that they are cocircular if they have the same 
center. A cocircular 4-component will be a 4-component with two cocircular opposite edges. 

Remark 3.7. It is easy to check that in a sequence of interior 4-components, if any of them is 
cocircular, then all the 4-components are. 

Lemma 3.8 ( W, Lemma 5.38]). Suppose we have a stationary graph with a sequence of at 
least three cocircular 4- components, so that the first and the last are boundary components and 
the remaining are interior ones {cocircularity refers to the edges of the boundary components 
meeting dD). Assume further that the components out of the chain belong to the same region. 
Then there is a continuous movement preserving perimeter and areas which creates an illegal 
meeting, so that the graph is not minimizing. 

Proof. Order the 4-components so that S7i and £l n are the boundary ones, and Qj meets fij+i. 
Let c, be the center of the cocircular arcs of flj. We can move these points, without changing 
neither d(ci, Cj+i) nor d(c±, 0) and d(c n , 0) in such a way that c\ and c n get closer and closer. 
With this movement of the points, we obtain a deformation of the graph which preserves 
perimeter and the areas of the regions and will create an irregular meeting. Hence the graph 
cannot be minimizing. □ 

Lemma 3.9 ( |W| Lemma 5.3]). Let e, f and g be three consecutive edges of a component, and 
let vi and v 2 be the corresponding vertices. Suppose e and g have the same geodesic curvature, 
and the angles in each vertex are the same. Let R be the line of points equidistant from v\ 
and V2- 

Then e and g are interchanged by the symmetry about R. Moreover, if e and g are cocircular, 
the common center lies in R, and if e and g are not cocircular, R coincides with the line of 
points equidistant from the centers of e and g. 

Given a graph C and a Jacobi function u defined on it, we will say that a point x in C is a 
nodal point if u(x) = 0. A nodal region will be a connected component of the complementary 
in C of the set of nodal points. 
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Proposition 3.10 ( HMRR Proposition 5.2]). Let C be a stationary graph separating the 
disk into three regions. Assume that there exists a Jacobi function with at least four nodal 
regions such that the nodal points are not vertices of the graph. Then C is unstable. 

Proof. Let u be the Jacobi function, and Ni, . . . ,N± nodal regions. Assume that the graph 
C is stable. For i = 1, . . . , 4 define Ui as the restrictions of u to Ni extended by zero to the 
whole graph. It is possible to find a nontrivial linear combination v of u±, U2, U3 so that the 
mean value conditions are satisfied for R±, R2 and R3. Moreover, v is an admissible 

function and, by stability, is a Jacobi function. As v vanishes on a subset of C containing N4 
and has nontrivial support contained in N± U N2 U N3, there is an edge t so that v vanishes on 
an open interval of £ but it is not identically zero on I. As v is a Jacobi function (a solution 
of a second order o.d.e.) this gives us a contradiction. □ 

Now we discard configuration [(8)| by a geometrical argument. The reader may compare 
this result with CHLL, Lemma 5.2.10] 

Proposition 3.11. Configuration\{8)\is not minimizing. 




Figure 5. Deformation of a chain of symmetric 4-components obtaining a 
graph with irregular vertices 

Proof. Suppose this configuration is minimizing. We will denote by C O2 C R2 the 
interior components, and by C R\, Q' 2 C R2 the boundary components. 

Consider any component of Ri or R2 and suppose the edges separating such a component 
and -R3 are cocircular. By Remark 13. 71 we will have a sequence of four cocircular 4-components 
and using Lemma 13.81 the configuration is not minimizing. 

Hence the considered edges cannot be cocircular, and applying Lemma l3.9| it is easy to 
check that there exists a horizontal symmetry of the chain of 4-components which is in fact 
a symmetry of the whole configuration. 

In each interior component, using Lemma 13.91 and taking into account that any pair of 
opposite edges cannot be cocircular, it can be seen that there exists a vertical symmetry, 
orthogonal to the horizontal one. 

If p\ = P2, we have that f2i and Q2 are congruent, so that they have the same area. 
Hence, we can interchange them, and after eliminating unnecessary edges, we obtain a new 
configuration enclosing the same areas with strictly less perimeter, which gives a contradiction. 
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If P2 = P3, applying Gauss-Bonnet Theorem to we obtain a contradiction. 
So we can suppose now p\ > P2 > P3- Let us distinguish two cases: 

Let z be the upper point in 80, 2 Pi 8D. Let 9 be the angle between the line Oz and a 
horizontal line. Assume that 9 is greater than or equal to tt/4. Then, the second coordinate 
of z will be greater than or equal to cos7r/4 = \[2j2. We can see 8Q, 2 as a vertical graph, 
so the length of 8Q' 2 will be greater than or equal to \[2. By similar arguments, taking into 
account that the point of 80,2 with maximum second coordinate will be higher than z, the 
length of 80,2 will be greater than or equal to 2^/2. Both quantities add up to more than 3, 
so by Lemma ll. 81 this configuration cannot be minimizing. 

Assume now that 9 is less than tt/4. Let I be the edge separating f2i and f2 2 , and v, v' its 
vertices. Consider a new edge I, the symmetric of / about the segment vv', and let T be the 
region enclosed between I and I. Move in the horizontal direction until one of its edges 
coincides with I. Reflect Q\ — T about a vertical axis so that the reflection of I coincides with 
the edge ^2 n (recall that each interior 4-component has a vertical symmetry, so the side 
edges have the same length and curvature). After eliminating an unnecessary edge, this new 
configuration will preserve length and areas, but it will be irregular. Hence, configuration [(8)| 
cannot be minimizing. It only remains to check that Q\ and wn l remain inside D under 
this geometrical transformation. 

In order to prove this it is enough to show that the portion of the original graph over the 
horizontal line L passing through z stays inside D. The portion £2 of the upper edge of 0,2 
over L makes an angle 9 with L, and can be translated horizontally to touch z at its boundary. 
Since 9 < tt/4, £2 must lie inside D. This implies that the transformation of O2 stays inside 
D. The upper edge £\ of Q\ has larger geodesic curvature than the upper edge of and 
makes a smaller angle with the horizontal line passing through v. This is enough to conclude 
that £\ must lie inside D. □ 



4. Proof of the theorem 

In this section we prove the main theorem. 

Theorem 4.1. Let C C D be a minimizing graph for three given areas. Then C is a unique 
standard graph. 

Proof. By Lemma 12.21 the graph C must be one of the listed in Figure 



Configurations |(1)| and (2) are unstable by Proposition 13.51 

Configuration [(3)| is also unstable: it is easy to check that the edges of Cyi are not cocircu- 
lar, so that by Remark 13 . 71 and Lemma 13.91 we have a vertical symmetry in this configuration. 
By Lemma 13.81 the edges of C13 cannot be cocircular if C is minimizing. Applying again 
Lemma 13.91 we get an horizontal symmetry. Both axes of symmetry will meet orthogonally at 
the origin, so we can consider the Killing field generated by rotations about 0, which yields 
a Jacobi function u vanishing on four points, one on the interior of each edge of the central 
4-component. So u has at least four nodal regions and we conclude by Proposition 13. 101 that 
C is unstable. 
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We now eliminate configurations |(4)| and |(5)| These configurations present an interior 
4-component of Ri, with three incident edges meeting the exterior of the disk. If we extend 
the fourth edge, it will meet dD orthogonally by Proposition an d we will obtain a config- 
uration of type ["(3) | Hence, the interior 4-component has two orthogonal symmetries meeting 
at the origin, and we conclude as before the existence of four nodal regions. 



Consider now configuration [(6)| Fix an interior 4-component O of R\. If we extend the 
edge leaving the boundary of Vt that does not reach dD, it will meet dD orthogonally due 
to the existence of a symmetry of which is in fact a symmetry of the disk. In this way we 
obtain a configuration of type |(3)] As above, Vl will have two orthogonal symmetries meeting 
at and so we can get four nodal regions yielding instability. 

Consider now configuration |(7)| Applying Lemma 1.3.91 the two 4-components will be 
symmetric about two lines r%, r<i passing through the center of the disk (the corresponding 
edges are not cocircular). Let q\, q2 be the intersection points of each line with the interior 
edges of these components, that will be zeros of the Jacobi field u induced by the one- 
parameter group of rotations about the origin. The reflection of q\ with respect to T2 lies in 
the boundary of the 3-component of -Ri and it is not a vertex of the configuration. This point 
is clearly also a zero of u. Then u has four nodal regions and the configuration is unstable. 

By Proposition 13 . 1 ll configuration [(8)| is not minimizing. 

Configuration [(9)| is unstable: if the top and bottom edges of each component are cocircular 
then the configuration is not minimizing by Lemma l3.8l Otherwise we can find an horizontal 
symmetry of the graph, which is also a symmetry of the disk by Lemma l3.9l Each interior 4- 
component has a vertical symmetry so that the interior components of Ri are identical. Using 
the function equal to +1 on one of these components, equal to —1 on the other component, 
and zero otherwise, we have obtained a function satisfying the mean value conditions (jl.lj) 
which is negative for the index form. So this configuration is unstable. We could also use the 



method of Proposition 13.111 to see that this configuration is nonminimizing. 

So the only remaining possibility is configuration [(10)1 the standard one. Uniqueness for 
given areas comes from Theorem 13.61 □ 



5. Final remarks 

In this paper we have obtained that the problem of dividing the disk into three areas has a 
unique solution in which all regions are connected, as in the problem of partitioning the disk 
into two areas. It is natural to conjecture that 

Conjecture 5.1. A minimizing graph separates the disk into connected regions. 

If we consider the problem for n regions, with n ^ 4, by Lemma l2.1l the region of largest 
pressure will have at most n — 1 nonhexagonal connected components and we can obtain by 
combinatorial arguments a list of all possible minimizing configurations. Of course the number 
of candidates increases very rapidly when the number of regions increases. We believe that 
the following conjectures are true 



Conjecture 5.2. The least perimeter way of dividing the unit disk into four regions of pres- 
cribed areas is given by configuration (1) of Figure El 
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Conjecture 5.3. The least perimeter way of dividing the unit disk into five regions of prescribed 
areas is given by configuration (2) of Figure El 




Figure 6. The conjectural configurations for n = 4 and n = 5 

For each case, we believe that there is another possibly stable configuration: for n = 4, the 
configuration with three boundary regions surrounding an interior one of three edges, and 
for n = 5, the one consisting in four boundary regions surrounding an interior region of four 
edges. But estimates we have done using Surface Evolver (Ken Brakke, 1992) for equal areas 
show that they are nonminimizing. Furthermore, for n = 4, if any of the areas tends to zero, 
we should obtain in the limit the standard configuration for three areas, which also discards 
the configuration described above at least for some areas. In the case n = 5 we should have 
the same behaviour. 

For n = 6 we give the following conjecture 

Conjecture 5.4. The least perimeter way of dividing the unit disk into six regions of prescribed 
areas is given by configuration of Figure [Tj 



Figure 7. The conjectural configuration for n = 6 



As before, we belive that the configurations of Figure |B1 below are stable, but estimates 
done with the Surface Evolver considering equal areas show that they are nonminimizing. 
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(1) (2) (3) 



Figure 8. Some other configurations for n = 6 
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